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ABSTRACT 

The distribution of 64 species of land snails (Mollusca: Gastropoda) across the 105 coun¬ 
ties of the State of Kansas was analyzed by principal components analysis to determine the 
relationship between species occurrence and climatic factors. Significant biotic responses of 
species assemblages to climatic gradients in precipitation and in annual temperature were dis¬ 
covered. In particular, total species diversity is positively correlated with precipitation (annu¬ 
al or summer alone) across the counties of Kansas, and the distributions of major clusters of 
gastropod species were correlated with average annual temperature (or latitude). 


Introduction 

This paper presents an analysis of the relationships 
between the distribution of species assemblages of land 
snails (Mollusca: Gastropoda) across the counties of the 
State of Kansas and gradients of climatic variables across 
those counties. For the large, mid-continental state of Kansas 
in the years represented by our data, the average annual pre¬ 
cipitation across the 105 geopolitical divisions (the counties) 
ranged from 15.82 to 41.68 inches, while the annual average 


temperature range was 50.6 to 58.9 °F, and the average num¬ 
ber of frost-free days varied between 154 and 198. 

Because the State of Kansas shows pronounced geo¬ 
graphic dines in climatic conditions, our aim was to explore 
whether there is a faunal response to these gradients in a 
group of organisms likely to be dependent upon climatic fac¬ 
tors. In a similar geographic study, Kadmon and Heller (1998) 
found a strong correlation between species distributions of 
land snails in Israel and the climatic gradients found there. 
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Materials and Methods 

The most recent published survey of species of land snails 
in Kansas remains the monograph of Leonard (1959), which 
gives distribution maps for each snail species in the counties 
of the state. Table 1 lists the 64 gastropod taxa found in 
Kansas, using nomenclature updated according to the system¬ 
atic treatment of Turgeon et al. (1988) and earlier works 
(Pilsbry, 1948; Taylor, 1960; Taylor, 1965). The distribution 
maps were converted into a presence/absence data table, 
where each of the 105 counties was scored according to 
whether each of the species was recorded as present. The 
counties of Kansas were coded according to a letter-and-num- 
ber checkerboard scheme (Figure 1), with lower letters of the 
alphabet to the east and higher numbers to the north. Thus, the 
most southeastern county (Cherokee) was coded Al, and the 
most northwestern county (Cheyenne) was coded 08. 

Climatic data for the sample localities were taken from the 
nearest reporting station of the United States Weather Bureau 
for the period 1899-1938, corresponding to the majority of 
the collection dates for the specimens collected by Leonard 
(1959) for the Museum of Natural History at the University of 
Kansas. The five climatic parameters that were used corre¬ 
spond to average values for the main weather station in each 
county, using published maps (Kincer, 1941). Their descrip¬ 
tions are given in Table 2, along with latitude and longitude. 

All statistical analyses were performed using the SPLUS 
2000 software package (MathSoft, 1999). We used principal 
components analysis (Morrison, 1967; Gauch, 1982) to reduce 
the complexity of the distributions of the 64 species across the 
105 counties in Kansas. This technique is commonly used in 
analyses of geographic distributions of biotic assemblages 
(Kadmon and Heller, 1998) to reduce the multivariate dimen¬ 
sionality of the original data. It constructs linear combinations 
of the original variables (in this case, the presences and 
absences of 64 species) that maximally encompass the varia¬ 
tion in the original data. In particular, the first principal com¬ 
ponent accounts for the greatest variance in the species data; 
the second component is perpendicular to and uncorrelated 
with the first component, and it accounts for the second largest 
amount of variation; the third component is orthogonal to the 
first two and accounts for the next largest amount of variation; 
etc. These synthetic variables thus represent uncorrelated 
aspects of species distributions across the counties of Kansas. 

The principal components for species distribution then were 
related to latitude, longitude, and the five specific climatic vari¬ 
ables for the 105 counties by using simple linear regression and 
correlation techniques (Gauch, 1982; Sokal & Rohlf, 1995). 

Results and Discussion 

Principal components analysis of distributions of 64 gas¬ 
tropod species 

Principal components analysis of the distributional data, 
based on the matrix of correlations among the 64 species 
variates, was successful in capturing a significant fraction of 



Figure 1. Map of the counties of Kansas, along with their let- 
ter-and-number codes, as follows: Al, Cherokee; A2, 
Crawford; A3, Bourbon; A4, Linn; A5, Miami; A6, 
Johnson; A7, Wyandotte; A8, Leavenworth; B1, Labette; B2, 
Neosho; B3, Allen; B4, Anderson; B5, Franklin; B6, 
Douglas; B7, Jefferson; B8, Atchison; B9, Doniphan; Cl, 
Montgomery; C2, Wilson; C3, Woodson; C4, Coffey; C5, 
Osage; C6, Shawnee; C7, Jackson; C8, Brown; Dl, 
Chautauqua; D2, Elk; D3, Greenwood; D4, Lyon; D5, 
Wabaunsee; D7, Pottawatomie; D8, Nemaha; El, Cowley; 
E2, Butler; E4, Chase; E5, Morris; E6, Geary; E7, Riley; E8, 
Marshall; F4, Marion; F5, Dickinson; F7, Clay; F8, 
Washington; Gl, Sumner; G2, Sedgwick; G3, Harvey; G4, 
McPherson; G5, Saline; G6, Ottawa; G7, Cloud; G8, 
Republic; HI, Haiper; H2, Kingman; H3, Reno; H4, Rice; 
H5, Ellsworth; H6, Lincoln; H7, Mitchell; H8, Jewell; II, 
Barber; 12, Pratt; 13, Stafford; 14, Barton; 16, Russell; 17, 
Osborne; 18, Smith; Jl, Comanche; J2, Kiowa; J3, Edwards; 
J4, Pawnee; J5, Rush; J6, Ellis; J7, Rooks; J8, Phillips; Kl, 
Clark; K2, Ford; K3, Hodgeman; K4, Ness; K6, Trego; K7, 
Graham; K8, Norton; LI, Meade; L2, Gray; L4, Lane; L6, 
Gove; L7, Sheridan; L8, Decatur; Ml, Seward; M2, Haskell; 
M3, Finney; M4, Scott; Nl, Stevens; N2, Grant; N3, Kearny; 
N4, Wichita; N6, Logan; N7, Thomas; N8, Rawlins; Ol, 
Morton; 02, Stanton; 03, Hamilton; 04, Greeley; 06, 
Wallace; 07, Sherman; 08, Cheyenne. 


the variation among all 105 counties in the first three prin¬ 
cipal components. Principal component 1, which accounted 
for 20.8 percent of the total variation in the data, loaded pos¬ 
itively on all 64 species, except Sitccinea pseudavcira, 
Succinea vaginacontorta , and Gastrocopta cristata; in other 
words, it contrasted the occurrences of these three species 
against those of the other 61 species, among which the high¬ 
est loadings were attributed to Zonitoides arboreus, 
Anguispira altemata, Philomycus carolinianus, Strobilops 
labyrinthicus, Gastrocopta armifera, Gastrocopta contrac- 
ta, and Carychium exile. The second principal component, 
accounting for another 11.2 percent of the variation, con¬ 
trasted 19 species against the other 45. The third component 
comprised an additional 6.2 percent of the variation and 
contrasted the distributions of a different set of 30 species 
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Table 1. Systematics of the gastropod species (or subspecies) found in Kansas. The 28 species subset used in a second analysis 
described in the text is indicated with boldface. 


Order Family Species 

Basommatophora 

Carychiidae 

Caryehimn exigmim 
(Say, 1822) 

Caryehimn exile exile 
I. Lea, 1842 

Stylommatophora 

Bulimulidae 

Biilimiilus dealbatus dealbatus 

(Say, 1821) 

Cochlicopidae 

Coclilicopa lubrica lubrica 
(Miller, 1774) 

Discidae 

Anguispira alternata altemata 
(Say, 1816) 

Haplotrematidae 

Haplotrema concavum 
(Say, 1821) 

Helicarionidae 

Euconidiis chersmus polygyrahis 
(Say, 1821) 

Helicodiscidae 

Heliocodiscus eigenmanni eigemnanni 
(Pilsbry, 1900) 

Heliocodiscus parallelus 
(Say, 1817) 

Heliocodiscus singleyanus singleyanus 
(Pilsbry, 1890) 

Philomycidae 

Philomycus carolinianus (Bose, 1802) 

Polygyridae 

Allogona profunda (Say, 1821) 
Euchemotrema leai aliciae 
(A. Binney, 1841) 

Mesodon clausus (Say, 1821) 
Mesodon inflectus inflectus (Say, 1821) 
Mesodon tliyroidus thyroidus 
(Say, 1816) 

Polygyra doifeailliana dorfeuilliana 
(I. Lea, 1838) 

Polygyra dorfeuilliana sampsoni 
(1. Lea, 1838) 

Polygyra jacksoni jacksoni 
(Bland, 1866) 

Stenotrema hirsutuin hirsutum 
(Say, 1817) 

Stenotrema stenotrema (Pfeiffer, 1842) 
Triodopsis albolabiis atteni (Say, 1816) 
Triodopsis cragini Call, 1886 
Triodopsis divesta (Gould, 1848) 
Triodopsis mnltilineata algonquinensis 
(Say, 1821) 

Triodopsis niultilineata multilineata 
(Say, 1821) 

Triodopsis neglecta (Pilsbry, 1899) 


Order 


Family 

Punctidae 

Pupillidae 


Succineidae 


Strobilopsidae 

Vitrinidae 

Zonitidae 


Valloniidae 


Species 


Punctual minutissimiun (1. Lea, 1841) 

Gastrocopta annifera (Say, 1821) 
Gastrocopta contracta (Say, 1822) 
Gastrocopta corticaria (Say, 1816) 
Gastrocopta cristata 
(Pilsbry & Vanatta, 1900) 
Gastrocopta holzingeri (Sterki, 1889) 
Gastrocopta pellucida liordeacella 
(Pfeiffer, 1841) 

Gastrocopta pentodon (Say, 1821) 
Gastrocopta procera procera 
(Gould, 1840) 

Gastrocopta tappaniana 
(C. B. Adams, 1842) 

Pupoides albUahris 
(C. B. Adams, 1841) 

Pupoides liordaceus (Gabb, 1866) 
Pupoides inomatus Vanatta, 1915 
Vertigo milium (Gould, 1840) 

Vertigo ovata ovata Say, 1822 
Vertigo tridentata Wolf, 1870 

Catinella vagans (Pilsbry, 1900) 
Catinella wandae (Webb, 1953) 
Oxyloma retuswn (1. Lea, 1834) 
Succinea concordialis Gould, 1848 
Succinea ovalis Say, 1817 
Succinea pseudavara Webb, 1954 
Succinea vagiuacontorta Lee, 1951 

Strobilops labyrintliicus (Say, 1817) 

Deroceras laeve (Muller, 1774) 

Hawaiia minuscula minuscuki 
(A. Binney, 1840) 

Mesompliinx enpreus ozarkensis 
(Rafmesque, 1831) 

Nesovitrea electiina (Gould, 1841) 
Nesovitrea identata identata (Say, 1823) 
Paravitrea capsella capsella 
(Gould, 1851) 

Paravitrea significans (Bland, 1866) 
Paravitrea simpsoni (Pilsbry, 1889) 
Striatura nieridionalis 
(Pilsbry & Ferriss, 1906) 

Striatura milium (E. S. Morse, 1859) 
Ventridens ligera (Say, 1821) 
Zonitoides arboreus (Say, 1816) 

Vallonia parvula Sterki, 1893 
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Figure 2. Location of the counties of Kansas (indicated with 
the locality codes of Figure 1) in the space defined by the 
first three principal components based on the distribution of 
64 gastropod species. 


against the remaining 34. Thus, the first three principal com¬ 
ponents together account for 38 percent of the variation in 
the distributional data of 64 gastropod species across the 
105 counties of Kansas. 

Inspection of the location of the individual counties in 
the three most important dimensions of principal compo¬ 
nents space (Figure 2) reveals that the most eastern coun¬ 
ties of Kansas (codes A1 to D8) occur at the high end of 
the first principal component, that the most southeastern 
county (Al: Cherokee) is distinct with respect to the sec¬ 
ond principal component, and that there is some separa¬ 
tion of the eastern counties, as well as H3 (Reno County) 
along the third component. The superposition of nearly all 
of the western counties indicate that they are distinct with 
respect to gastropod assemblages. 


Table 2. The seven climatic variables and correlations of 
the first three principal components (comps.) of species 
distributions of 64 gastropod species with these climatic 
variables. The last five variables represent annual averages. 



Comp. 1 

Comp. 2 

Comp. 3 

Latitude (°N) 

0.151 

0.303 

-0.115 

Longitude (°W) 

-0.730 

-0.040 

0.099 

Annual precipitation (inches) 

0.730 

0.010 

-0.132 

Summer precipitation (inches) 

0.721 

0.034 

-0.150 

Annual temperature (°F) 

0.125 

-0.203 

-0.023 

July temperature CF) 

-0.128 

-0.091 

-0.040 

Number of frost-free days 

0.436 

0.157 

0.013 


Table 3. Correlations between principal component 
(comp.) scores based on 28 gastropod species and climatic 
variables in the counties of Kansas. 


Species 

Comp. 1 

Comp. 2 

Comp. 3 

Euchemotrema leai aliciae 

0.241 

-0.136 

0.001 

Mesodon thyroidus thyroidus 

0.139 

-0.009 

-0.332 

Mesodon clausus 

0.163 

0.202 

0.092 

Triodopsis albolabris alleni 

0.219 

0.213 

-0.073 

Bnlimulus dealbatus dealbatus 

0.170 

-0.122 

-0.104 

Euconulus chersinus polygyratus 

0.230 

0.011 

-0.097 

Nesovitrea electrina 

0.244 

0.136 

-0.048 

Nesovitrea identata identata 

0.249 

-0.033 

-0.085 

Paravitrea significcins 

0.082 

-0.030 

-0.321 

Hawaiia minuscnla minuscula 

0.179 

-0.180 

0.166 

Striatura milium 

0.091 

0.385 

0.351 

Striatum meridionalis 

0.091 

0.386 

0.350 

Zonitoides arboreus 

0.247 

-0.108 

-0.011 

Anguispira altemata altemata 

0.258 

0.061 

-0.130 

Heliocodiscus parallelus 

0.180 

-0.236 

0.149 

Punctum minutissimnm 

0.167 

0.238 

0.053 

Philomycus carolinianus 

0.204 

0.028 

-0.266 

Succinea concordialis 

0.110 

0.141 

0.161 

Catinella vagans 

0.070 

0.186 

0.169 

Strobilops labyrinthicus 

0.254 

0.102 

-0.195 

Gastrocopta armifera 

0.188 

-0.254 

0.199 

Gastrocopta contracta 

0.244 

-0.113 

-0.009 

Gastrocopta pentodon 

0.189 

0.099 

-0.174 

Gastrocopta tappaniana 

0.160 

-0.089 

0.066 

Gastrocopta procera procera 

0.151 

-0.222 

0.262 

Pupoides albilabris 

0.167 

-0.230 

0.199 

ValIonia parvula 

0.157 

-0.282 

0.286 

Carychium exile exile 

0.214 

0.247 

-0.013 


Relationships between species distributions and climatic 
gradients 

Table 2 gives the correlations between the first three prin¬ 
cipal components and the values for each of the seven cli¬ 
matic factors for each county. It is clear, for example, that the 
first component is strongly positively correlated with total 
precipitation (r = 0.730, P < 0.001) and with summer pre- 
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cipitation (r = 0.721, P < 0.001), but strongly negatively cor¬ 
related with longitude (r = 0.730, P < 0.001). As might be 
expected from the mid-continental location of Kansas, there 
is a strong positive correlation between annual and summer 
precipitation (r = 0.997, P < 0.001). Furthermore, this state 
partly was chosen for this study because of the strong rela¬ 
tionship (r = 0.981, P < 0.001) between longitude and pre¬ 
cipitation, with a remarkably linear gradient from the wetter 
east to the drier west (see Figure 3). However, Figure 4 
demonstrates that annual precipitation is a good predictor of 
the first principal component only for the eastern counties, 
where there generally is a higher diversity of gastropod 
species (Figure 5). A much better predictor of the first com¬ 
ponent (r = 0.976, P < 0.001) is simply the species diversity 
(Figure 6), although as noted above not all species load pos¬ 
itively on this component. Furthermore, the supeiposition of 
many of the western counties along the first component that 
was evident before can be explained, in large part, by their 
depauperate gastropod populations. 

The second principal component is most strongly posi¬ 
tively correlated with latitude (r = 0.303, P < 0.01) and neg¬ 
atively correlated (r = 0.203, P < 0.01) with annual temper¬ 
ature (Table 2). Since higher latitudes generally correspond 
with lower temperatures over the counties of Kansas (r = 
0.800, P < 0.001), these relationships imply that the second 
most important axis of gastropod species distribution can be 
explained partially by a second, less pronounced climatic 
gradient. The third principal component correlates only 
weakly with climatic variables; the highest correlation is 
with summer precipitation, but it is not statistically signifi¬ 
cant (r = 0.150, P > 0.05). 

Principal components analysis of distributions of 28 
selected gastropod species 

In order to test the robustness of these results, a second prin¬ 
cipal component analysis was undertaken using a subset con¬ 
sisting of 28 of the 64 gastropod species found in Kansas. Table 
1 indicates in boldface type the species selected for this analy¬ 
sis, which represent most of the families of gastropods. They 
include most of the species that loaded highly in the analysis of 
the full dataset, as listed above. Figure 7 shows the disposition 
of the 105 counties in the space defined by the first two princi¬ 
pal components. Again, a longitudinal trend is evident (with 
most of the western counties superimposed), plus the separa¬ 
tion of B4 (Anderson County) along the second component. 

The first principal component in this case loaded positively 
on all 28 species (see Table 3), with the highest contributions 
provided by Angitispira ciltematci, Strobilops labyrinthicus , 
Nesovitrea indentata, Zonitoides arboreus , Nesovitrea electri- 
na, Gastrocopta contracta, and Euchemotrema leai. Many of 
these same species contributed highly to the first component of 
the full 64-species dataset. Once again, scores on the first com¬ 
ponent across the counties of Kansas were highly positively 
correlated with the climatic factors of annual precipitation (r - 
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Figure 3. Relationship between longitude and total annual 
precipitation in Kansas. 



Figure 4. Relationship between total annual precipitation 
and scores on the first principal component of gastropod 
species distributions in Kansas. 



Figure 5. Relationship between longitude and gastropod 
species diversity in Kansas. 
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Figure 6. Relationship between gastropod species diversity 
and the first principal component. 



Figure 7. Location of the counties of Kansas in the space 
defined by the first two principal components based on the 
distribution of 28 gastropod species. 
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Figure 8. Relationship between total annual precipitation 
and scores on the first principal component of 28 gastropod 
species distributions in Kansas. 


0.784, P < 0.001; Figure 8) and summer precipitation (r = 
0.781, P < 0.001), were highly negatively correlated with lon¬ 
gitude (r = 0.773, P < 0.001), and were very strongly correlat¬ 
ed with species diversity ( r = 0.966, P < 0.001). 

The second principal component most strongly contrasted 
the distributions of Striatum meridionalis. Striatum milium, 
Carychium exile , and Punctum minutissimum with those of 
Vallonia parvula, Gastrocopta armifera, Helicodiscus paral- 
lelus , and Pupoides albilabris. This component was significant¬ 
ly correlated only with latitude (r = 0.238, P < 0.05). The third 
component contrasted the distributions of Striatum meridion¬ 
alis, Striatum milium, Vallonia parvula, and Gastrocopta pro- 
cera with those of Mesodon thyroidus, Paravitrea significans, 
and Philomycus carolinianus. It was correlated significantly 
only with two climatic variables: positively with latitude (r = 
0.372, P < 0.01) and negatively with average temperature (r = 
0.231, P < 0.05); of course, these two variables are correlated 
negatively with each other, as noted above. 

The congruence of the results for the full set of 64 gastropod 
species found in Kansas and the selected subset of 28 species 
strengthens our interpretation that there has been a significant 
biotic response of this group of gastropods to climatic gradients 
in precipitation (from east to west across the checkerboard of 
Kansan counties) and, less strongly, in annual temperature (from 
south to north). These conclusions are comparable to the simi¬ 
larly strong response found by Kadmon and Heller (1998) of the 
land snails of Israel to a strong rainfall gradient from north to 
south. To our knowledge, the present study represents the first 
analysis of the relationship between land snail distributions and 
climatic factors in the central region of the United States. 

Acknowledgements 

We are grateful for the helpful comments of J. Heller 
(Hebrew University of Jerusalem) and J. B. Keiper (The 
Cleveland Museum of Natural History) on the manuscript. 
We also appreciate the financial support of Cleveland State 
University for publication costs, through a supplemental 
EFFRD grant to M. Tevesz. 

References 

Gauch, H. G., Jr. 1982. Multivariate Analysis in Community 
Ecology. Ordination, Chapter 4, p. 109-172. Cambridge 
University Press, Cambridge, England. 

Kadmon, R., and J. Heller. 1998. Modelling faunal responses to cli¬ 
matic gradients with GIS: land snails as a case study. Journal of 
Biogeography, 25:527-539. 

Kincer, J. B. 1941. Climate and weather data for the U.S., p. 685-741. 
In G. Hambridge (ed.). Climate and Man (USDA Yearbook), U.S. 
Government Printing Office, Washington. 

Leonard, A. B. 1959. Handbook of Gastropods of Kansas. University 
of Kansas Museum of Natural History, Miscellaneous Publication 
20. Lawrence. 225 p. 

MathSoft 1999. SPLUS 2000. Data Analysis Products Division, 
MathSoft, Inc., Seattle. 








2002 


CLIMATIC GRADIENTS 


11 


Morrison, D. F. 1967. Multivariate Statistical Methods. The 
structure of multivariate observations: I. Principal compo¬ 
nents, Chapter 7, p. 221-258. McGraw-Hill, New York. 

Pilsbry, H. A. 1948. Land Mollusca of North America. American 
Academy of Natural Sciences, Monograph 3, 2:521-1113. 

Sokal, R. R., and F. J. Rohlf. 1995. Biometry: The Principles and 
Practice of Statistics in Biological Research. Linear regres¬ 
sion, Chapter 14, p. 451-544, and Correlation, Chapter 15, p. 
545-608. W. H. Freeman & Co., New York. 

Taylor, D. W. 1960. Late Cenozoic molluscan faunas from the 
High Plains. United States Geological Survey Professional 
Paper 337. 94 p. 


Taylor, D. W. 1965. The study of Pleistocene nonmarine mollusks 
in North America, p. 597-611. In H. E. Wright and D. G. Frey 
(eds.), The Quaternary of the United States. Princeton 
University Press, Princeton. 

Turgeon, D. D., A. E. Bogan, E. V. Coan, W. K. Emerson, W. G. 
Lyons, W. L. Pratt, C. F. E. Roper, A. Scheltema, F. G. 
Thompson, and J. D. Williams. 1988. Common and Scientific 
Names of Aquatic Invertebrates from the United States and 
Canada: Mollusks. American Fisheries Society Special 
Publication 16. American Fisheries Society, Bethesda. 277 p. 


